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Abstract 

We have extended our program of QCD simulations with an improved Kogut-Susskind quark 
action to a smaller lattice spacing, approximately 0.09 fm. Also, the simulations with a « 0.12 
fm have been extended to smaller quark masses. In this paper we describe the new simulations 
and computations of the static quark potential and light hadron spectrum. These results give 
information about the remaining dependences on the lattice spacing. We examine the dependence 
of computed quantities on the spatial size of the lattice, on the numerical precision in the com- 
putations, and on the step size used in the numerical integrations. We examine the effects of 
autocorrelations in "simulation time" on the potential and spectrum. We see effects of decays, 
or coupling to two-meson states, in the ++ , 1 + , and 0~ meson propagators, and we make a 
preliminary mass computation for a radially excited 0~ meson. 

PACS numbers: 11.15Ha,12.38.Gc 
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I. INTRODUCTION 



We have extended our ongoing program of lattice QCD simulations with three flavors 
of dynamical quarks. In this paper we describe the new simulations we have done, and 
present spectrum results for the light hadrons and the static quark potential. In a previous 
workjl| we presented results for these quantities from a set of runs with a lattice spacing 
of approximately 0.12 fm and light quark masses ranging down to 0.2 times the estimated 
strange quark mass. Since that time we have extended the a ~ 0.12 fm runs to smaller 
quark masses, and increased the statistics on the m u ^ = 0.2 m s run. More importantly, we 
have done simulations at a smaller lattice spacing of approximately 0.09 fm in quenched 
QCD and with three dynamical flavors at three values of the light quark mass: m u ^ = m s , 
mu,d = 0.4m s and m u ^ = 0.2 m s where m s is the strange quark mass estimated before doing 
the simulations 2]. This enables us to address the question of lattice spacing effects, i.e., 
extrapolation to the continuum, to greater accuracy than we could before. Two short runs 
were made at larger integration step size than used in the main simulation as an additional 
check on the systematic errors in the simulation algorithm. At our smallest quark mass, we 
have computed the hadron propagators in double precision on a subset of the lattices as a 
check on the numerical accuracy of the computations. Finally, we have done an explicit test 
of the effects of the finite spatial size of the simulated system by adding a run with a larger 
spatial size than in the main run. 

In addition to the light hadron spectrum, the gluon configurations generated in this 
program are being used for computations of the static quark potential [3|], heavy quark 
and heavy-light meson spectroscopyO, [J , heavy-light meson deca y co nstants Q, f n , fx, 
and chiral 0(p 4 ) parameters j^f Q, |J a s ^, exotic meson masses icl . the topological sus- 
ceptibility in QCd[tl1]. masses [3, |13|, |1J|, and parton 
distributions [1^]. For those quantities where accurate lattice results are available and sys- 
tematic errors are relatively well understood, there is good agreement with experimental 
values among a large set of quantities [3]. While this work focuses on describing the simula- 
tions, the static potential, and the light hadron spectrum, results from these other quantities 
are important in our analysis. In particular, the T mass splittings give the most accurate 
estimates of the lattice spacing, and several of these quantities enter into our estimates of 
the correct strange quark mass. In turn, some of the results presented here, such as the 
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dependence of the static potential on the lattice spacing, and the tests of the effects of 
molecular dynamics step size and spatial size of the lattices, are important in evaluating 
these other works. 



II. SIMULATIONS 

The simulations used here are a continuation those described in Ref. which contains 
a more detailed description of the simulation program. We use an improved Kogut-Susskind 
quark action, the "a 2 ^' or "Asqtad" action, which removes lattice artifacts up to order a 2 g 2 . 
Configurations were generated using the hybrid-molecular dynamics "R algorithm" ^|, with 
separate pseudofermion fields for the light and strange quarks, except where all three quarks 
are degenerate. The momenta conjugate to the gauge fields were refreshed at the end of 
every trajectory, with the trajectory length being one simulation time unit. Lattices were 
archived every six time units, and the hadron spectrum and static quark potential were 
calculated on these stored lattices. 

Table U] summarizes the parameters of the runs. For completeness, it includes runs re- 
ported in Ref. , although we will not repeat tabulation of masses from runs that have not 
been extended since that time. In identifying runs, we will quote the light (degenerate u 
and d) and strange quark masses as ami/ s = 0.01/0.05, for example. 

III. STATIC POTENTIAL AND LENGTH SCALE 

We use the static quark potential to relate the lattice spacings in our different runs. In 
particular, we use the quantity r\ defined by r\F{rx) = 1.00. We choose r x because of 
its ease and accuracy of computation and lack of dependence on the valence quark mass. 
Computation of this quantity and the effects of dynamical quarks on the potential have been 
discussed in Refs. 0, [J . Here we add points at smaller quark mass and, more importantly, 
points at a finer lattice spacing which allow a preliminary continuum extrapolation. As 
before, we fit to the form in Ref. [lsj ]. 



V{r) — C + ar — a/r + A (Vf ree (r) - 1/r) . (1) 

where Vf ree (r) is the potential calculated in free field theory, using the improved gauge action. 
This lattice correction term is used at distances less than 3a. 
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am Ujd 1 am s 


10/ ff 2 


L 


u 


res. 


e 


lats. 


a/n 


quenched 


8.00 


20 


0.8879 


na 


na 


408 


0.3762(8) 


0.02 / na 


7.20 


20 


0.8755 


1 x 10" 


-4 


0.013 


370 


0.3745(14) 


0.40 / 0.40 


7.35 


20 


0.8822 


2 x 10" 


-5 


0.03 


332 


0.3766(10) 


0.20 / 0.20 


7.15 


20 


0.8787 


5 x 10" 


-5 


0.03 


341 


0.3707(10) 


0.10 / 0.10 


6.96 


20 


0.8739 


5 x 10" 


-5 


0.03 


339 


0.3730(14) 


0.05 / 0.05 


6.85 


20 


0.8707 


1 x 10" 


-4 


0.02 


425 


0.3742(15) 


0.04 / 0.05 


6.83 


20 


0.8702 


5 x 10" 


-5 


0.02 


351 


0.3765(14) 


0.03 / 0.05 


6.81 


20 


0.8696 


5 x 10" 


-5 


0.02 


564 


0.3775(12) 


0.02 / 0.05 


6.79 


20 


0.8688 


1 x 10" 


-4 


0.0133 


484 


0.3775(12) 


* 0.01 / 0.05 


6.76 


20 


0.8677 


1 x 10" 


-4 


0.00667 


658 


0.3852(14) 


* 0.01 / 0.05 


6.76 


28 


0.8677 


1 x 10" 


-4 


0.00667 


241 


0.3814(14) 


* 0.007 / 0.05 


6.76 


20 


0.8678 


1 x 10" 


-4 


0.005 


493 


0.3783(13) 


* 0.005 / 0.05 


6.76 


24 


0.8678 


5 x 10" 


-5 


0.003 


197 


0.3796(19) 


* quenched 


8.40 


28 


0.8974 


na 




na 


396 


0.2681(5) 


* 0.031 / 0.031 


7.18 


28 


0.8808 


2 x 10" 


-5 


0.02 


496 


0.2613(9) 


* 0.0124 / 0.031 


7.11 


28 


0.8788 


5 x 10" 


-5 


0.008 


527 


0.2698(9) 


* 0.0062 / 0.031 


7.09 


28 


0.8782 


5 x 10" 


-5 


0.004 


592 


0.2714(9) 



TABLE I: Parameters of the improved action simulations. A "*" at the beginning of the line 
indicates a run which is new or has been extended since the report in Ref lj. The first column 
gives the light and strange quark masses in lattice units, and the second column, the gauge coupling. 
"L" is the spatial size of the lattice. The time size is 64 for the coarse lattices and 96 for the fine 
lattices, uq is obtained from the average plaquette. The conjugate gradient residual tabulated here 
is the residual used in generating configurations; a smaller residual was used in computing hadron 
propagators, "e" is the time step size in configuration generation. The second to the last column 
is the number of stored lattices, and the last column is the lattice spacing in units of n determined 
from the static potential in this run. A "smoothed" lattice spacing, discussed later, will be used 
to convert results to physical units. The last four lines, with a ~ 0.09 fm, will be referred to as 
"fine" lattices. 
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FIG. 1: A "shape parameter" for the static potential, tq/t\. The red octagons are from coarse 
(a 0.12 fm) lattices with three degenerate quark flavors, and the red squares from coarse lattices 
with two light and one strange quark. At (M^/Mp) 2 = 0.15 the upper square is from the L = 28 
run and the lower from the L = 20 run. The blue crosses are from the fine (a ~ 0.09 fm) runs. The 
single green diamond is from a two flavor simulation. The magenta burst is the continuum and 
chiral extrapolation discussed in the text, with the smaller error bar the statistical error and the 
larger the systematic error. In this figure we have chosen to use (M^/Mp) 2 for the abscissa instead 
of the (Mn-ri) 2 used in other figures because this lets us put the entire range of quark masses up 
to the quenched limit (M^ — > oo) in the graph. 
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am u>d 1 am s 


10/5 2 


r\ja (run) 


r\ja (smoothed) 


0.0492/0.082 


6.503 


1.774(10) 


1.778 


0.0328/0.082 


6.485 


1.786(10) 


1.788 


0.0164/0.082 


6.467 


1.783(12) 


1.797 


0.0082/0.082 


6.458 


1.807(10) 


1.802 


0.082/0.082 


6.561 


1.816(10) 


1.805 


0.0492/0.0492 


6.475 


1.807(28) 


1.766 


0.0328/0.0328 


6.470 


1.768(30) 


1.828 


0.0164/0.0164 


6.430 


1.796(22) 


1.813 


0.0492/0.0492 


6.500 


1.818(23) 


1.821 


0.0492/0.0492 


6.450 


1.735(30) 


1.713 


0.0328/0.0328 


6.450 


1.757(30) 


1.784 


0.0164/0.0164 


6.450 


1.857(25) 


1.858 


0.0082/0.0082 


6.420 


1.843(20) 


1.827 


0.005/0.050 


6.76 


2.634(13) 


2.632 


0.007/0.050 


6.76 


2.644(09) 


2.623 


0.010/0.050 


6.76 


2.598(08) 


2.610 


0.010/0.050 


6.76 


2.621(09) 


2.610 


0.020/0.050 


6.79 


2.649(08) 


2.650 


0.030/0.050 


6.81 


2.656(10) 


2.662 


0.040/0.050 


6.83 


2.666(11) 


2.673 


0.050/0.050 


6.85 


2.679(11) 


2.683 


0.030/0.030 


6.79 


2.678(14) 


2.650 


0.031/0.031 


7.18 


3.827(12) 


3.822 


0.0124/0.031 


7.11 


3.707(13) 


3.711 


0.0062/0.031 


7.09 


3.687(12) 


3.684 



TABLE II: Smoothed r±/a compared with r±/a determined from each run. The top block is from 
lattices with a « 0.18 fm from tuning runs for our high temperature simulations, while the second 
and third blocks are the "coarse" and "fine" lattices respectively. Five short "tuning runs" are 
omitted from this table. Several of the runs have been extended since fitting of the smoothed r\ 
was done. 



While we expect r\ /a to be a smooth function of the quark masses and gauge couplings, 
r\ja determined from fitting the potential in a particular run will have a statistical error, and 
fluctuate from its ideal (infinite statistics) value. To minimize the effects of these run-to-run 
fluctuations, we have fit a smoothed ri/a for our three flavor lattices with quark masses less 
than or equal to the strange quark mass. Over the range of masses and gauge couplings we 
have used, a simple fitting form 

login/a) = C 00 + C w - 7.0^ + On (2m M + m s ) + C 20 - 7.0^ (2) 

gives an acceptable fit with a x 2 of 30.3 with 26 degrees of freedom, with 

C 00 = 1.2578(27) 

C w = 0.9371(93) 

C 01 = -0.828(29) 

C 20 = -0.271(22) (3) 



Table Uj] shows values of r\/a used in the fit together with the smoothed r\/a for each run. 
We have used this smoothed r\ja in converting results from units of the lattice spacing into 
units of T\. 

The shape of the static quark potential is affected by dynamical quarks. One of many 
possible ratios parameterizing this shape is the ratio r /r 1 . We use the results in Fig. [TJto 
extrapolate ro/r\ to the physical quark mass and continuum limit. Simultaneously fitting 
coarse and fine lattice results to a constant plus linear terms in the quark mass and a 2 a s 
gives 

tq/tx = 1.476(7) -0.049(10)(M^/M P ) 2 -0.12(4)(a/ri) 2 a s (a)/a s (0.12fm) , (4) 

with x 2 = 3-6 for 8 degrees of freedom, using a s from Ref. jjj. In fitting the potential the 
same distance range, y/2 — 6, was used for all the coarse lattices, and range Vo — 7 for all 
the fine lattices. Therefore, the statistical error bars in Table |H] and Fig. ^ appropriately 
represent the fluctuations in r\fa or ro/rx within each of these two sets of runs. However, 
there is a systematic effect from the choice of fit range which is common to all coarse runs and 
all fine runs, but may differ between the two sets. Varying the fitting range over reasonable 



8 



ranges suggests that this systematic error can be conservatively estimated as an uncertainty 
of 0.01 in the difference between the coarse and fine lattice Tq/tx. This leads to a systematic 
uncertainty of about 0.018 in the continuum extrapolation, leading to an estimate 

ro/rx = 1.474(7)(18) (5) 

at the physical M n /M p in the continuum limit. 

To compute t\ in physical units, we need to set the lattice scale using a directly measurable 
physical quantity. A convenient choice is the T spectrum, in particular the 2S-1S and 1P- 
1S splittings. This gives a scale a -1 = 1.588(19) GeV on the coarse 0.01/0.05 lattices, and 
a- 1 = 2.271(28) GeV on the fine 0.0062/0.031 lattices Q- For light quark masses <m s /2, 
the mass dependence of these quantities and of r\ appears to be slight, and we neglect it. 
With our smoothed values of r^/a, we then get r\ = 0.324(4) fm on the coarse lattices and 
ri = 0.320(4) fm on the fine lattices. 

To extrapolate r% to the continuum, we first assume that the dominant discretization 
errors go like a^a 2 . Using oty(q*) jj], |3] (with scale q* = 3.33/a) for as gives a ratio 
(a sa 2 )fi ne / \oisa 2 ) coarse = 0.428. Extrapolating away the discretization errors linearly then 
results in r\ = 0.317(7) fm in the continuum. However, taste- violating effects, while formally 
0(a 2 s a 2 ) and hence subleading, are known to be at least as important as the leading errors 
in some cases. Therefore, one should check if the result changes when the errors are assumed 
to go like a 2 s a 2 . Taking as = «y(3.33/a) gives a ratio (a|a 2 )fi ne /(a|a 2 ) coarse = 0.375; while 
a direct lattice measurement of the taste-splittings to be presented in the next section gives 
a ratio of 0.35. Extrapolating linearly to the continuum then implies T\ = 0.318(7) fm or 
ri = 0.319(6) fm respectively, in agreement with the previous result. For our final result, 
we use an "average" ratio of 0.4 and add the effect of varying this ratio in quadrature with 
the statistical error. We obtain n = 0.317(7) (3) fm. The second error is a crude estimate 
of the systematic error from the choice of fit ranges for the static potential. 

A similar calculation to estimate r yields 0.471(6) fm on the coarse run and 0.466(6) fm 
on the fine run, with a continuum extrapolated value of 0.462(H)(4) fm, where the second 
error is an estimate of the systematic error from choice of fit ranges in the potential. If we 
take the above estimate of tq/ti and multiply by r\ = 0.317 fm, we obtain instead ro = 0.467 
fm, and the difference in these two calculations of tq is another measure of systematic error. 
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IV. LIGHT HADRON MASSES 

Our procedures for calculating and fitting hadron propagators are described in Ref. 
With the exception of the non-Goldstone pions at am u ^ = 0.0124, we used Coulomb gauge 
wall sources, with eight source time slices evenly spread through the lattice. Propagators 
were fit with varying minimum distances, and with the maximum distance either at the mid- 
point of the lattice or where the fractional statistical errors exceeded 30% for two successive 
time slices. In most cases, to reduce the effect of autocorrelations, propagators from four 
successive lattices (24 simulation time units) were blocked together before computing the 
covariance matrix. Masses were selected by looking for a combination of a "plateau" in the 
mass as a function of minimum distance and a good confidence level (x 2 ) for the fit. We also 
made an effort to choose minimum distances that are smooth functions of the couplings, 
recognizing that statistically we should have some fits with low and high confidence levels. 

A. Pseudoscalar mesons 

We calculated masses for the exact Goldstone (75 ® 75) pseudoscalar mesons in all of 
the runs. For the ami/ s = 0.0124/0.031 run we calculated the masses of all of the different 
taste pions, allowing us to see how the taste symmetry breaking decreases with lattice size. 
Figure El shows the fitted masses for the pion, the kaon and the "unmixed ss" from the fine 
lattice run with m Uj d = 0.2 m s . Table ITTT1 shows the selected fits for the pseudoscalar meson 
masses. 

With Kogut-Susskind quarks there are four "tastes" of valence quark, and hence sixteen 
different tastes of pseudoscalar mesons, grouped in eight multiplets. In the continuum limit 
these are degenerate, and the improved action reduces these splittings relative to the one-link 
fermion action. In our previous work on the coarse lattices we verified that these pion masses 
show the partial taste symmetry restoration predicted by Lee and Sharpej^j]. In particular, 
we expect near degeneracy between pairs of pions between which 70 is replaced by 7^, e.g. 
taste 7075 with taste 7^75. Also, the squared masses are approximately linear in the quark 
mass, with all tastes having the same slope. This means that a dimensionless measure of 
taste symmetry breaking, (M 2 — Mq) r 2 , is almost independent of the quark mass. Having 
verified these properties on the coarse lattice, we computed non-pointlike pion propagators 
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Q-TTlvalence 




amps 


range X^/D conf. 


0.015 (vr) 


OO 


0.21643(14) 


18-47 25/28 0.62 


0.03 (vr) 


oo 


0.30259(14) 


24-47 21/22 0.53 


0.01 (vr) 


0.01/0.05 


0.22439(20) 


19-31 9.1/11 0.61 


0.01 (vr) 


0.01/0.05 


0.22421(12) 


19-31 4.7/11 0.94 


0.007 (vr) 


0.007/0.05 


0.18881(19) 


20-31 14/10 0.18 


0.005 (vr) 


0.005/0.05 


0.15970(20) 


22-31 11/8 0.19 


0.01/0.05 (K) 


0.01/0.05 


0.38327(22) 


17-32 23/14 0.067 


0.01/0.05 (K) 


0.01/0.05 


0.38304(20) 


17-32 14/13 0.38 


0.007/0.05 (K) 


0.007/0.05 


0.37268(25) 


20-31 8.6/10 0.57 


0.005/0.05 (K) 


0.005/0.05 


0.36550(29) 


20-31 6.4/10 0.78 


0.05 (ss) 


0.01/0.05 


0.49427(18) 


17-32 19/14 0.18 


0.05 (ss) 


0.01/0.05 


0.49443(18) 


17-31 17/13 0.20 


0.05 (ss) 


0.007/0.05 


0.49317(19) 


20-31 12/10 0.31 


0.05 (ss) 


0.005/0.05 


0.49276(23) 


20-31 5.2/10 0.87 


0.031 (vr) 


0.031/0.031 


0.32003(18) 


25-47 20/21 0.52 


0.0124 (vr) 


0.0124/0.031 0.20638(18) 


30-47 22/16 0.15 


0.0062 (vr) 


0.0062/0.031 0.14794(19) 


35-47 7/11 0.8 


0.0124/0.031 (K) 0.0124/0.031 0.27209(18) 


30-47 23/16 0.11 


0.0062/0.031 (K) 0.0062/0.031 0.25319(19) 


30-47 14/16 0.61 


0.031 (ss) 


0.0124/0.031 0.32585(17) 


27-47 29/19 0.07 


0.031 (ss) 


0.0062/0.031 0.32727(14) 


32-47 5.6/14 0.97 



TABLE III: Pseudoscalar meson masses. Here we include runs that are new or have been extended 
since Ref. jjj; results at larger quark masses can be found there. The first column is the valence 
quark mass(es), and the second column the sea quark mass or masses. The particle name is in 
the first column. Here "vr" indicates valence quark mass equal to the lighter dynamical quarks, or 
degenerate in the quenched case. U K" indicates one valence quark equal to the light dynamical 
quarks and one at about m s , while "ss" indicates a fictitious meson with two valence quarks with 
mass about m s , in a flavor nonsinglet state. The remaining columns are the hadron mass, the 
time range for the chosen fit, x 2 an d number of degrees of freedom for the fit, and the confidence 
level of the fit. The first block is from the quenched run at W/g 2 = 8.4, the second block from 
the coarse three flavor runs, and the last block from the fine three flavor runs. The two lines with 
am sea = 0.01/0.05 are from the runs with L = 20 and 28. 
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FIG. 2: Pseudoscalar masses as a function of minimum distance included in the fit from the 
run with 10/<? 2 = 7.09 and amij s = 0.0062/0.031. The size of the symbols is proportional to the 
confidence level of the fit, with the size of the symbols in the labels corresponding to 50%. These 
fits included only a single exponential. Fits selected to quote in the mass tables are marked with 
arrows. 

on only one of the fine lattice runs, with lO/V = 7.11 and 

wmi/s = 0.0124/0.031, which 

has a lattice spacing of a/r x = 0.269. In Table HV1 we give these pion masses, together with 
those from the coarse lattice run with comparable quark masses. To facilitate comparison, 
these masses are given in units of r±. We also give the measure of taste symmetry breaking, 
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(M 2 — Mq) r\, for these masses. It can be seen that (M 2 — M G ) r\ for each taste on the fine 
lattices is consistently about 0.35 times the value on the coarse lattices. This is consistent 
with the expected scaling as a 2 a 2 s described above, which, using a s = ay(q*) and q* = 3.33/a 
suggests a ratio of 0.375. 



pion taste 


M n r\ (coarse) 


(Ml - Mr) r? (coarse) 


M n r% (fine) 


(Ml - Ml) r\ (fine) 


ratio 


/ 


0.8251(45) 




0.7659(7) 






7075 


0.9386(19) 


0.2003(35) 


0.8127(11) 


0.0739(18) 


0.369(11) 


7i75 


0.9426(16) 


0.2078(30) 


0.8116(26) 


0.0721(42) 


0.347(21) 


7i7j 


1.0033(34) 


0.3259(69) 


0.8372(41) 


0.1143(68) 


0.351(22) 


7i7o 


1.0044(29) 


0.3280(59) 


0.8383(26) 


0.1162(44) 


0.354(15) 


7i 


1.0555(53) 


0.4334(12) 


0.8576(56) 


0.1489(95) 


0.344(22) 


7o 


1.0558(32) 


0.4339(67) 


0.8602(37) 


0.1534(64) 


0.354(16) 


1 


1.1029(80) 


0.5358(75) 


0.8899(93) 


0.2054(165) 


0.383(31) 



TABLE IV: Taste symmetry violations on coarse and fine lattices. The second and fourth columns 
contain the masses for the different pions in units of r\ for a coarse and fine lattice run. The coarse 
lattice run (from Ref. ]|) was at 10/g 2 = 6.79 and ami/ s = 0.02/0.05, and had a lattice spacing 
a/ri = 0.377. The fine lattice run was at W/g 2 = 7.11 and am^ s = 0.0124/0.031, and had a 
lattice spacing ajr\ = 0.269. The physical quark masses are similar, as evidenced by the similar 
Goldstone pion masses. The third and fifth columns are a measure of taste symmetry breaking, 
(M 2 — Mq) rl, on the coarse and fine lattices, and the final column is the ratio of this measure 
between the fine and coarse lattice runs. 



In a separate analysis we calculate "partially quenched" pseudoscalar masses and decay 
constants, where the valence quark and sea quarks have different masses j?], 0. These re- 
sults have been analyzed using chiral perturbation theory including terms parameterizing 
the taste symmetry breaking|22|. From this analysis we find and fx at the physical 
quark masses, and values for several of the low energy constants in chiral perturbation 
theory. Another product of the computations of nip 5 and fps is a determination of the 
lattice quark masses corresponding to the real world. We define the strange and light quark 
masses at fixed lattice spacing, am l f" and am^, to be the lattice masses that give the ex- 
perimental values for M K and M n . To determine am[ at and am l ^ d , we fit the mass and 
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decay constant data to chiral log forms that take into account staggered taste violations 



2J. We find am l s at = 0.0390(1) (1^), am l ^ d = 0.00148(1) (tt) on the coarse lattices, and 
am iat = o.0272(1)(+Jq), am x ^ d = 0.00103(0)(4) on the fine lattices, where the errors are sta- 
tistical and systematic. The systematic error is dominated by that coming from the chiral 
extrapolation/interpolation and the ~2% scale uncertainty. 

We have also calculated masses of excited pseudoscalar mesons. Because this requires 
consideration of two-meson states, discussion of this is deferred to a later section on hadronic 
decays and excited states. 



B. Vector mesons 



Figure El shows vector meson masses versus minimum distance fit for the fine lattice run 
with the lightest quark mass. Mass estimates for all of the runs are in Table Note 
that despite our relatively small quark masses, none of these vector mesons are below the 
threshold for decay into two pseudoscalars, since the angular momentum of the vector mesons 
requires that the vector meson at rest decay into pseudoscalars with momentum 2tt/L. In 
addition, we require a combination of tastes in the pseudoscalars that overlaps with the 
taste of the vector meson — the vector mesons tabulated here have spin <g> taste = 7« ® 7i- 

Table IV II shows masses for 1 + mesons. These mesons can decay into a vector and a 
pseudoscalar meson, and these simulations reach into the quark mass region where this 
threshold is crossed. We defer discussion of this effect to the next section. 



C. Baryons 

Table IVTT1 contains masses for the octet nucleon and 5. We do not tabulate the A and £ 
since our code does not cleanly separate the light quark isospins. In principle, the nucleon 
mass could be fit by methods similar to those used for the pion mass and decay constant, 
incorporating effects of continuum chiral corrections, lattice artifacts like taste symmetry 
breaking, finite size effects and partial quenching. Such an analysis is not yet available. 
However, statistical errors on the nucleon mass are much larger than for the pseudoscalars, 
so this full machinery may be less important here. An alternative strategy for dealing 
with lattice artifacts is to perform a continuum extrapolation at the quark masses used 
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Q-TTlvalence 




amy 


range 


X 2 /D 


conf. 


0.015 (p) 


oo 


0.4660(30) 


11-25 4/11 


0.97 


0.03 (p) 


oo 


0.4992(15) 


5-25 


18/15 


0.28 


0.01 (ft) 


0.01/0.05 


0.5690(50) 


6-22 


15/13 


0.32 


0.01 {p) 


0.01/0.05 


0.5680(30) 


6-19 


10/10 


0.42 


0.007 (p) 


0.007/0.05 


0.5510(40) 


6-18 


11/9 


0.26 


0.005 (p) 


0.005/0.05 


0.5340(80) 


6-15 


5.8/6 


0.44 


0.01/0.05 (#*) 


0.01/0.05 


0.6492(25) 


8-23 


5.2/12 0.95 


0.01/0.05 (K*) 


0.01/0.05 


0.6462(18) 


8-27 


29/16 


0.023 


0.007/0.05 (K*) 


0.007/0.05 


0.6330(30) 


9-23 


10/11 


0.54 


0.005/0.05 (K*) 


0.005/0.05 


0.6190(40) 


10-23 


15/10 


0.12 


0.05 (6) 


0.01/0.05 


0.7193(14) 


9-30 


11/18 


0.90 


0.05 (0) 


0.01/0.05 


0.7194(11) 


9-31 


15/19 


0.74 


0.05 (0) 


0.007/0.05 


0.7114(16) 


12-30 12/15 


0.69 


0.05 (4>) 


0.005/0.05 


0.7140(30) 


14-29 15/12 


0.25 


0.031 (p) 


0.031/0.031 


0.4781(14) 


16-42 36/23 


0.043 


0.0124 (p) 


0.0124/0.031 0.4173(13) 


10-33 31/20 


0.059 


0.0062 (p) 


0.0062/0.031 0.3895(28) 


10-27 11/14 


0.65 


0.0124/0.031 (K*) 


0.0124/0.031 0.4483(18) 


15-42 42/24 


0.013 


0.0062/0.031 (K*) 


0.0062/0.031 0.4350(11) 


10-34 13/21 


0.91 


0.031 ((f>) 


0.0124/0.031 0.4831(8) 


14-47 55/30 


0.0032 


0.031 (0) 


0.0062/0.031 0.4810(40) 25-45 18/17 


0.39 



TABLE V: Vector meson masses. Runs tabulated and the format are the same as in Table Mil 
Here "p" indicates valence quark mass equal to the lighter dynamical quarks, or degenerate in the 
quenched case, "if*" indicates one valence quark equal to the light dynamical quarks and one 
at about m s , while "0" indicates two valence quarks with mass about m s , although in a flavor 
nonsinglet state. 
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FIG. 3: Vector meson masses as a function of minimum distance included in the fit from the run 
with W/g 2 = 7.09 and am l/s = 0.0062/0.031. 

in simulations, and then fit these extrapolated masses to continuum chiral perturbation 
theory. Figure H] shows the nucleon masses in units of r x . This graph also contains a very 
rough sketch of how such a continuum and chiral extrapolation might begin. The right 
most magenta fancy plus is a linear extrapolation in a?a s of the coarse and fine results at 
m u,d ~ 0.4m s to a = 0, as indicated by the red line. The middle fancy plus is a similar 
continuum extrapolation at m u ^ ~ 0.2 m s . The solid straight line is a linear extrapolation 
to the physical pion mass. As a rough estimate of the effects of chiral logarithms, the two 



16 



am* 



ampy range X^/D conf. 



0.015 (01) oo 
0.03 (a x ) oo 
0.015 (b 
0.03 (fei) 



1 oo 



oo 



0.720(40) 9-25 11/11 0.48 

0.730(6) 7-25 10/13 0.67 

0.741(22) 6-25 7.3/14 0.92 

0.748(10) 7-25 15/13 0.33 



0.01 (ai) 0.01/0.05 

0.01 (en) 0.01/0.05 

0.007 (ai) 0.007/0.05 

0.005 (oi) 0.005/0.05 

0.01 (6i) 0.01/0.05 

0.01 0.01/0.05 

0.007 0.007/0.05 

0.005 (h) 0.005/0.05 



0.820(40) 6-15 5.3/6 0.50 

0.848(24) 6-17 6.4/8 0.60 

0.767(21) 5-15 10/7 0.16 

0.790(40) 5-15 9.6/7 0.21 

1.020(90) 6-22 15/13 0.32 

0.980(60) 6-19 10/10 0.42 

0.810(40) 5-18 11/10 0.34 

0.700(90) 6-15 5.8/6 0.44 



0.031 (ai) 0.031/0.031 0.667(4) 8-25 11/12 0.56 

0.0124 ( ai ) 0.0124/0.031 0.600(8) 8-30 22/19 0.30 

0.0062 (ai) 0.0062/0.031 0.532(19) 10-26 14/13 0.36 

0.031 (bx) 0.031/0.031 0.681(5) 7-25 21/13 0.08 

0.0124 (6i) 0.0124/0.031 0.632(9) 7-33 34/23 0.07 

0.0062 (6i) 0.0062/0.031 0.650(50) 10-27 11/14 0.65 



TABLE VI: Pseudovector meson masses. Runs tabulated and the format are the same as in 
Table ITTT1 



curved lines are chiral perturbation theory forms constrained to match the two continuum 
extrapolated points. These forms have two free parameters, so we emphasize that this is 
not a fit and there is no test of consistency of these forms with our data. The yellow upper 
curved line is an expansion in powers of M n up to order M% log(M vr ) from Ref. and the 
cyan lower curve is a form where the nucleon-delta mass splitting is also treated as small j^. 
It is clear that fine lattice results at a smaller quark mass will be needed, since the slopes of 
the chiral perturbation theory forms are clearly different from the lattice results for quark 
masses as small as 0Am s . 
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Q-TTlvalence 




arris 


range X 2 /D conf. 


0.015 (N) 


oc 


0.6267(18) 


8-23 14/12 0.30 


0.03 (N) 


oo 


0.7134(18) 


12-30 11/15 0.78 


0.01 (N) 


0.01/0.05 


0.7710(40) 


6-16 4.4/7 0.74 


0.01 (N) 


0.01/0.05 


0.7670(30) 


6-17 3.8/8 0.88 


0.007 (N) 


0.007/0.05 


0.7480(30) 


5-14 5/6 0.54 


0.005 (N) 


0.005/0.05 


0.7230(60) 


5-14 9.8/6 0.13 


0.01/0.05 (S) 


0.01/0.05 


0.9810(30) 


8-20 5.3/9 0.81 


0.01/0.05 (S) 


0.01/0.05 


0.9737(20) 


8-21 16/10 0.09 


0.007/0.05 (3) 


0.007/0.05 


0.9670(50) 


10-20 9.2/7 0.24 


0.005/0.05 (5) 


0.005/0.05 


0.9540(80) 


10-19 4.6/6 0.60 


0.031 (N) 


0.031/0.031 


0.6996(11) 


7-37 50/25 0.0023 


0.0124 (N) 


0.0124/0.031 0.5815(19) 


10-29 17/16 0.41 


0.0062 (N) 


0.0062/0.031 0.5190(40) 


11-23 4.6/9 0.87 


0.0124/0.031 (3 


) 0.0124/0.031 0.6696(17) 


13-33 12/17 0.80 


0.0062/0.031 (3 


) 0.0062/0.031 0.6519(18) 


12-30 19/15 0.21 



TABLE VII: Octet baryon masses. Runs tabulated and the format are the same as in Table IrTTI 

V. TESTS OF SYSTEMATIC AND STATISTICAL ERRORS 

The results in the previous two sections allow us to make several algorithm tests as well 
as more physical tests. 

A. Single versus double precision 

As the valence quark masses are made smaller, the condition number of the fermion 
matrix increases and one might worry that double precision is necessary for computing 
the hadron propagators. In general, we have used single precision for the computations 
at each lattice site, with global sums in double precision. At our smallest quark mass, 
o^u.d — 0.005, we have tested the accuracy of our hadron spectrum and static potential 
computations by repeating the computation in double precision on a subset of the lattices. 
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0.0 0.5 1.0 1.5 2.0 

(m„ r t ) 2 

FIG. 4: Nucleon masses. The blue diamonds and octagons are quenched coarse and fine runs 
respectively. The red squares are three flavor coarse lattice results, and the red bursts the three 
flavor fine lattices. The magenta fancy plusses connected by the straight line and the two curved 
lines are continuum and chiral extrapolations discussed in the text. The fancy diamond is the 
experimental value, with an error bar from the uncertainty in r\. 

Table IVlTTI shows results for a number of quantities evaluated on a set of 137 lattices with 
ami/ s = 0.005/0.050. Note that since these are fit on exactly the same sets of lattices 
with exactly the same programs, any discrepancies are the result of the different precision. 
However, we provide statistical errors to show how the effects of roundoff compare with the 
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statistical errors. For all of these quantities the effects of using single precision are small 
compared with the statistical errors, and with the statistical errors we would get from any 
reasonable lengthening of this run. 



Quantity 


Double 


Single 


Comment 


V(2,0,0) 


0.829883(852) 


0.829888(853) 


potential at r=(2,0,0) 

JT \ 7 7 / 


F(2,2,2) 


1.05426(503) 


1.05451(502) 




F(3,3,3) 


1.2511(194) 


1.2511(194) 




ri/a 


2.63933(1679) 


2.63915(1678) 


t=4-5, block=5 


U / -L 


1.4566(64) 


1.4566(64) 


t=4-5, block=5 


vr(20) 


411.53(1.55) 


411.44(1.55) 


prop, at d=20 


p(6) 


143.76(1.78) 


143.73(1.78) 




aM n 


0.15965(22) 


0.15966(21) 


d=20-31, x 2 /D = 0.60 


aM K 


0.36519(34) 


0.36519(34) 


d=20-32, X 2 /D = 0.82 


aM p 


0.5330(83) 


0.5330(83) 


d=6-14, X 2 /D = 0.85 


aM N 


0.7311(84) 


0.7312(84) 


d=6-14, X 2 /D = 0.50 



TABLE VIII: Comparison of results with single and double precision computations. The first three 
lines are the static quark potential at three different spatial separations. These separations are in 
the spatial region used in fitting the potential. The next two lines are parameters extracted from 
fitting the potential, the inverse lattice spacing in units of n and a "shape parameter" ro/n. The 
second part of the table contains hadron propagator comparisons. The 7r(20) and p(6) show the 
pion and rho propagators summed over a time slice at time separations 20 and 6. These distances 
are near the minimum of the range used in fitting the masses, and so are among the most important 
distances in our fits. Finally, the last four lines are hadron masses computed from the double and 
single precision propagators. 

B. Integration step size 

Our simulation algorithm is expected to introduce errors proportional to e 2 where e is 
the simulation time step size. Based on previous experience and our expectations about the 
scaling of the fermion force with the quark mass, we have used a step size of about 2/3 of 
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the light quark mass in these runs. As a check on these effects, we have made short runs 
with larger step sizes at one of our small quark masses (the same parameters at which we 
checked effects of the spatial size of the lattice.) The production runs here were done at a 
step size of e = 0.0067 (658 lattices), and the short tests at step sizes of 0.01 (49 lattices) 
and 0.01333 (53 lattices) with lattice size 20 3 x 64. Table 11X1 shows results for the static 
quark potential and some hadron masses at these different step sizes, using the same fitting 
ranges in each case. Since the short runs were too short for a good error analysis, statistical 
errors on these quantities are estimated by scaling the errors on the L = 20, e = 0.0067 run 
by the square root of the ratio of the numbers of configurations used. 





L = 20 


L = 28 


L = 20 


L = 20 


Q. 


e = 0.0067 


e = 0.0067 


e = 0.0100 


e = 0.013 


□ 


1.70092(2) 


1.70094(3) 


1.70096(7) 


1.70066(7) 




0.07421(10) 


0.07420(13) 


0.07374(37) 


0.07488(35) 


n/a 


2.598(8) 


2.621(9) 


2.649(29) 


2.619(28) 


aM w 


0.22439(20) 


0.22421(12) 


0.22500(73) 


0.22554(70) 


aM p 


0.569(5) 


0.568(3) 


0.557(18) 


0.558(18) 


aM N 


0.771(4) 


0.767(3) 


0.785(15) 


0.753(14) 



TABLE IX: Effect of integration step size. These are from runs with W/g 2 = 6.76 and am^ s = 
0.01/0.05. Columns two and three are our long runs with L = 20 and 28 using a step size of 0.0067. 
(Our usual practice is to use a step size about 2/3 of the lightest quark mass.) Columns four and 
five are from short runs with step sizes 0.01 and 0.01333. 

C. Spatial size of the lattice 

In one of our coarse lattice runs, 10/g 2 = 6.76, ami/ s = 0.01/0.05, we have made a second 
run at a larger spatial lattice size, 28 3 x 64. (We have also lengthened the run with L = 20, 
so this is the run where we have the best statistics.) This allows us to explicitly check the 
effects of the spatial lattice size. Table IXl shows the results of this test for the static quark 
potential and simple hadron propagators. Note that these values of r 1 /a fall on opposite 
sides of the interpolated ("smoothed ri") value of 2.610, and the values of r /ri fall on 

21 



opposite sides of a straight line fit to the coarse lattice points in Fig. [TJ leading us to believe 
that we do not see any statistically significant finite size effects in either the potential or 
the hadron masses. The sizes of these two lattices in physical units are 2.43 and 3.40 fm, 
using r\ = 0.317 fm to set the scale, and M^L is 4.48 and 6.27 respectively. Using the 
(staggered) chiral fits [?], 0] to light pseudoscalar masses and decay constants, it is possible 
to estimate the leading finite volume correction on M n . We expect a difference A = 0.00026 
between L = 20 and L = 28 results, consistent with the observed value in the simulations, 
A = 0.00018(23), shown in Table H 



Quantity 


L = 20 


L = 28 


A 


n/a 


2.598(8) 


2.621(9) 


-0.023(12) 


ro/n 


1.4461(36) 


1.4533(34) 


-0.0072(50) 


M w 


0.22439(20) 


0.22421(12) 


0.00018(23) 


M p 


0.569(5) 


0.568(3) 


0.001(6) 


M N 


0.771(4) 


0.767(3) 


0.004(5) 



TABLE X: Comparison of results with different spatial sizes. These are from the runs with W/g 2 = 
6.76 and om^ s = 0.01/0.05. The spatial sizes were L = 20 and 28, corresponding to physical 
sizes of 2.4 and 3.4 fm, using n = 0.317 fm to set the physical scale. The first two lines are 
parameters extracted from fitting the potential, the inverse lattice spacing in units of r\ and a 
"shape parameter" ro/n. The second part of the table contains hadron mass comparisons. A is 
the L = 20 value minus the L = 28 value in each row. 



D. Autocorrelations 

Because of the high cost of generating sample configurations with dynamical quarks, suc- 
cessive samples were taken at simulation time intervals such that they are not completely 
statistically independent. The resulting autocorrelations (in simulation time) affect the sta- 
tistical errors on all of the computed quantities. The "exponential autocorrelation time", 
which is determined by the eigenvalue of the Markov process matrix which is closest to one, is 
expected to be the same for all calculated quantities. However, the contribution of this slow- 
est mode to various quantities varies, and to parameterize the effect of autocorrelations on 
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individual quantities we use the "integrated autocorrelation time", r int = J2 s Cq(s), where 
s runs over the simulation time separations and Cq(s) is the normalized autocorrelation for 
quantity Q, 

r ( , _ m + s)Q(t))-(Q)i 

Cq[ ) ~ mmrw ■ () 

Because we need a covariance matrix to calculate masses from the average propagators, and 
getting a nonsingular covariance matrix requires more samples than there are points in the fit 
range, we cannot get a hadron mass from one sample. So, to study autocorrelations of hadron 
mass estimates we use the "mirror image" of this procedure — we do single elimination 
jackknife fits with one sample omitted from the data set and compute the autocorrelations of 
these jackknife fits. Figure shows the jackknife pion masses as a function of the simulation 
time of the omitted sample for the run with 10/g 2 = 6.76 and ami/ s = 0.01/0.05. For 
example, Table IXll shows Cq(6) where Q is the tt, p or nucleon mass or the amplitude in the 
pion propagator, and the simulation time separation is six units, corresponding to successive 
stored lattices. From this table we can see that the normalized autocorrelation is largest 
for the pion mass, and has no obvious systematic dependence on the light quark mass. 
Therefore, we average the autocorrelations over the quark masses, separately for the coarse 
and fine runs. The resulting autocorrelations as a function of simulation time separation are 
plotted in Fig. El 

Not surprisingly, the autocorrelation times are larger on the fine lattices than on the 
coarse lattices. In Ref. ^j] autocorrelations of the topological charge were computed on these 
lattices. The topological charge evolves more slowly than the hadron masses, with estimated 
autocorrelation times as lar ge a s 35 time units for the 10/g 2 = 7.18, am^ = 0.031/0.031 



run. We refer the reader to 



111 ] for more discussion. 



VI. HADRONIC DECAYS AND EXCITED STATES 

When the quark mass is small enough, most of the hadrons we study are unstable, 
decaying strongly into two or more lighter hadrons. In principle, although not always in 
practice, fitting to the ground state mass in our propagators will give the mass of the 
lightest state with the right quantum numbers in the periodic box, which in many cases will 
be a two particle state. In Ref. Q we showed this effect in the ++ (do) channel. Figure [3 
updates this plot with more results on coarse lattices at light quark mass, and the new 
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io/s 2 




N 


M n 


An 


M p 


M N 


6.85 


0.05 


425 


0.196 


0.079 


0.047 


0.077 


6.83 


0.04/0.05 


351 


0.383 


0.127 


-0.031 


0.119 


6.91 


0.03/0.05 


564 


0.274 


0.161 


0.082 


0.070 


6.79 


0.02/0.05 


486 


0.173 


0.169 


0.025 


0.143 


6.76 


0.01/0.05 


658 


0.229 


0.056 


0.046 


0.014 


6.76 


0.007/0.05 


487 


0.150 


0.056 


-0.055 


-0.020 




average 


2971 


0.229 


0.106 


0.024 


0.062 


7.18 


0.031 


496 


0.426 


0.223 


0.074 


0.203 


7.11 


0.0124/0.031 


534 


0.311 


0.142 


-0.002 


0.034 


7.09 


0.0062/0.031 


586 


0.283 


0.152 


0.055 


0.011 




average 


1616 


0.336 


0.170 


0.042 


0.078 



TABLE XI: Normalized autocorrelations Cq(6) for hadron masses and the pion amplitude in the 
light quark runs. The third column is the number of samples in each run. We also show the results 
averaged over all the coarse runs, and over all the fine runs, where the third column is the total 
number of coarse or fine lattices. 

results on the fine lattices. For the three flavor runs, the fine lattice points agree well with 
the coarse lattice results. The figure also shows the mass of the lowest energy two-meson 
state expected to couple to this particle, n + rj. Surprisingly, the new points at the lighest 
quark masses increasingly deviate from this two-meson mass, which is not understood. The 
light mass quenched propagators remain difficult to fit, which may not be surprising for 
unstable particles in an unphysical theory. We have not yet tried fitting to the particle-plus- 
ghost form suggested by Bardeen et al. j^. For quark masses where the two-meson state 
has lower energy, it would be satisfying to find a one meson (do) state as an excited state in 
the propagator. Our attempts to do this have been unsuccessful so far. In the fine lattice 
run at ami/ s = 0.0062/0.031 we were able to extract an excited state mass, shown as the 
cyan decorated square in Fig. [7| However, the mass of this state is still much smaller than 
the extrapolations from large quark mass, and it is likely also a two-meson state, perhaps 
KK. 

We also expect to see the pseudovector (1 + ) mesons couple to two zero-momentum 
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FIG. 5: Single elimination jackknife masses for the pion, from the run with IO/3 2 = 6.76 and 
am z/s = 0.01/0.05, using fits with D m \ n = 19. 

mesons, although for these mesons we are not as far below the threshold as in the ++ 
case. Figure |H1 shows 1 H (61) masses as a function of quark mass along with the decay 
channel mass M p + M^. We tentatively attribute the downturn at the lightest quark masses 
to this decay, although better statistics at the lightest coarse lattice and a lighter mass fine 
lattice run would clarify the situation. Again, we are unable to get good fits for the lightest 
mass quenched propagators. 

Kogut-Susskind meson propagators generally include normal exponential contributions 
from one J PC value and an oscillating exponential component from a parity partner state. 
In the case of the Goldstone pion, the parity partner has the exotic J PC = + ~ and thus 
does not contribute to the propagator. In combination with a relatively high signal-to-noise 
ratio at all time separations, this enhances our ability to determine the 0~ + contributions. 
Specifically, in addition to the one-state fits, which we presented in Figure El and Table ITTT| 
when we performed a two-state fit of the pseudoscalar propagator data, we were able to 
determine the mass of a second, excited 0~ + state. We have presented preliminary results 
of this analysis in 2(|. We fit _+ propagators to the form: 



C(t) = A {e- Mot + e 



-Mo(T-t) 



)+A 1 (e- Mlt + e 



-Mi(T-t) 



(7) 
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FIG. 6: Normalized autocorrelations for the w, p and nucleon masses and the amplitude of the 
pion propagator as functions of the separation in simulation time. Results for all the coarse lattice 
runs are averaged together, as are all the fine lattice runs. 

where A and M are the amplitude and mass of the ground state, and A\ and M\ are 
the same for the lowest excited state. Figure El is a sample pion fit plot showing the fitted 
values of clMq and aM\ as a function of the minimum time separation, D m [ n , included 
in the fits. By comparing to one-state fits shown in Figure El note the inclusion of an 
excited state in the fitting function allows high-confidence fits to extend down to a D min 
of 2 or 3, as might be expected. The excited state's contribution to the propagator decays 
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FIG. 7: ++ energies. The red squares are three flavor coarse runs, and the red fancy plusses, three 
flavor fine runs. The blue octagons are a quenched coarse run and the blue burst a quenched fine 
run. The cyan decorated square is an excited ++ mass from one of the runs. The green diamonds 
and crosses are sums of ir and r\ masses on coarse and fine runs respectively, where the r\ mass is 
estimated from M 2 = + §m 2 j, with m S s the unmixed ss pseudoscalar mass. The straight 
blue line is an extrapolation of ao masses from heavier quark runs (not shown in this graph). 

to unresolvable levels relatively quickly, however, and consequently larger fit distances are 
often not so useful. Figure El summarizes the two-state fits for the _+ masses as a 
function of (M^ri) 2 . These excited state masses fit a linear function of (M^ri) 2 to a 12% 
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FIG. 8: 1^ energies. The red squares are three flavor coarse runs and the red bursts three 
flavor fine runs. The blue octagons and crosses are quenched coarse and fine runs respectively. 
The green diamonds and fancy plusses are sums of rho and pion masses on coarse and fine three 
flavor runs respectively. The fancy diamond on the left is the experimental value, with an error 
bar corresponding to the uncertainty in r\. 

confidence level. As the statistical errors on the excited pion mass fits are large compared 
with the differences between the coarse and fine lattice fits, we considered all of the mass fits 
together in the linear fit. Extrapolating the resulting linear function to the physical value of 
(M^n) 2 = 0.050, we get a prediction of a physical 0~+ excited state at 1362(41) (247)MeV, 
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FIG. 9: Fit plot showing ground state and excited state masses of Goldstone pion ^ as a 
function of -Dminj the minimum distance included in the fit. This is from the run with W/g 2 = 7.09, 
ami/s = 0.0062/0.031, with fits using D max = 28. The symbol size is proportional to confidence 
level. 



which agrees within the large errors with the mass of the 7r(1300) state. The first error is 
statistical. The second is the systematic error predominantly due to contributions to the 
propagator which are unaccounted for in the form of the fitting function. We estimate this 
by examining the fit plots and estimating the range of mass values one might reasonably 
choose, that is, this error reflects the stability of the fitted value under variation of the fit 
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FIG. 10: Ground state and excited state + masses as a function of (M^ri) 2 . In the legend, the 
symbol on the left represents fine lattice results, and that on the right the coarse lattice results. 

range, e.g., the difference between the D min = 3 and D min = 4 points in Figure M and is 
reflected in Figures El and E] as light cyan error bars on the excited states. We linearly 
extrapolate the individual systematic errors to (M n ri) 2 = 0.05. Systematic errors due to 
chiral extrapolation, finite lattice size and lattice spacing, are small relative to the statistical 
error and the systematic error from additional states. 

Similarly, an excited state is evident in the h ss propagator. The analysis of states con- 
taining strange quarks is complicated by the fact that our simulated strange quark masses, 
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am s = 0.050, 0.031 differ from the physical strange quark mass, am* = 0.039, 0.027 (for the 
coarse and fine lattices respectively) as discussed in subsection IIV Al To correct for this, 
after fitting to the form of Eq. Q , we interpolated the meson masses to the correct physical 
values of the strange quark mass, m*, using 

Mp S (m s ) = Mp S (m s ) - (m s - m s ) , (8) 

m s - m u4 

where we use the mass of the excited 0~ + state at the simulation value of m s for Mps(m s ), 
and the pion excited state on the same lattices for M-p${m u ^) ■ We cannot interpolate masses 
from lattices with three flavors of degenerate quarks in this manner, so we eliminate them 
from this analysis. 

The interpolated excited state masses fit a linear function of (M n ri) 2 and we again ex- 
trapolated the resulting form to the physical (M n ri) 2 . The result is M S g = 1645(40) (145) 
for the excited s~s psuedoscalar state. 

We have no pure ss physical 0~ + with which to compare ground state fits. We can, how- 
ever, compare the extrapolation of the corrected excited state masses with the experimental 
mass of the ?y(1440), which one expects to be dominated by the ss contribution. This is 
consistent with our result with the large systematic error. We display all of the pion and 
(corrected) _+ s~s fits in Figure E3 with physical states for comparison. 

Even more interesting is the kaon propagator. Formed of a light quark and a strange 
quark, the kaon, J p = 0~, has no definite charge-conjugation quantum number when m u ^ ^ 
m s . Consequently, it has a non-exotic parity partner with J p = + , and the propagator 
has a tiny, but significant oscillating component. On theses lattices the amplitude of the 
oscillating state is significantly smaller than that of the kaon ground state, and the mass 
is greater than that of the kaon ground state, thus it does not interfere with with single- 
exponential fits of the propagator at large time separations (D m m > 14). Two-state fits to 
the form of Eq. (J7J) fail at all time separations because the + mass falls below that of the 
first excited 0~ state. Figure fTTI shows an attempt to fit the 10 / g 2 = 7.09, am u ^ = 0.0062, 
am s = 0.031 fine lattice propagator to two non-oscillating exponentials, as in Eq. (J7J). All 
fits are of extremely low confidence levels and there is no evident plateau for the excited 
state. Figure shows fits of the same propagator to a three-state form, 

C(t) = A (e- Mot + e~ M ^) + A!(e- Ml * + e"^^) + A 2 {-l)\e- M * + e~ M ^), (9) 
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with high confidence levels and masses of consistent value through a large variation in the 
lower limit of the fit range, -D m in- 

Propagators from both fine lattice sets with m u ^ 7^ m s were inconsistent with double 
exponential forms, (Eq. |7j), but fit to triple exponentials, (Eq. |UJ) with high confidence. The 
same was true of the coarse lattice sets with am u ^ < 0.02. In general, we find that as 
^u,d — ► 777-5, the amplitude of the oscillating state becomes indistinguishable from zero, 
presumably because charge conjugation regains its status as a good quantum number. In 
the fits to kaon propagators from the coarse lattice set with 10/g 2 = 6.79, am u ^ = 0.030 
we were no longer able to distinguish the amplitude of any oscillating state from noise. 
Confidence levels for both two-state and three-state fits were a few tenths of a percent, yet 
we could discern equivalent plateaus for the excited 0~ state mass as a function of D min in 
each case. Attempts to read a plateau for the oscillating state were unconvincing. For the 
coarse lattices with 10/ g 2 = 6.81, am u ^ = 0.040, and both coarse and fine lattices with three 
degenerate flavors of quarks, two state fits resolved the excited state with high confidence 
(as we have mentioned before when we considered these very same fits as limiting cases of 
both pions and 0~ + ss mesons.) 

The oscillating scalar state is far lighter than the lightest strange + meson, the Kq(1A30). 
It does, however, agree well with the sum of the masses of the dominant _K"q(1430) decay 
mode products, K + tt, on every lattice set for which it was measured. Resolution of the 
Kq(1430) — > K + tt decay channel is additional evidence that our simulations with light 
dynamical quarks correctly reproduce the expected complexities of the physical world. When 
we perform similar fits to quenched kaon propagators we can find no evidence of an oscillating 
+ state, even with widely separated valence quark masses, such as ami/ s = 0.0062/0.031. 
Furthermore, with the quenched kaon propagators, it is simple to extract the contribution 
of the first excited 0~ state, see for example Fig. [T3J 

We have also performed an extrapolation of the excited kaon state masses to the physical 
value of (M^ri) 2 = 0.050. Again considering the fine and coarse lattice data together the 
excited states fit, with 8% confidence level, to a line which intercepts (M n ri) 2 = 0.050 at 
1527(46) (68)MeV. This is in decent agreement with the K(1460) state and inconsistent with 
the K(1460)'s expected decay products, ttttK, which should be at about 775 MeV. This 
lends credence to the belief that the K(1460) is a true mesonic state. 

Figure El summarizes the fits to the kaon propagators. As with the ss states, we have 
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corrected the ground state and excited state mass fits for the difference between the simulated 
strange quark mass and the physical strange quark mass using the interpolation expression 
(JSJ). Since we have measured a + state at only one value of the strange quark mass for each 
lattice spacing, interpolation of the + state is not possible. We include the pion ground state 
and the sum of the pion and (uncorrected) kaon ground state masses for comparison. We 
include isospin-averaged physical states for comparison. We display these results numerically 
in Table IXTTl 



W/g 2 


am U)d /am s 


•^states 


aM + 


^0+ 


aM n + aM K 




range 


conf 


6.85 


0.05 


2 






0.97 


1.05(2)(10) 


3-18 


0.36 


6.83 


0.04/0.05 


2 






0.90 


1.02(3)(2) 


4-32 


0.36 


6.81 


0.03/0.05 


2 






0.81 


1.07(3)(5) 


4-26 


0.008 


6.79 


0.02/0.05 


3 


0.63(12)(10) 


-3(2) 


0.72 


0.96(3)(2) 


3-16 


0.39 


6.76 


0.01/0.05 


3 


0.76(15)(4) 


-13(9) 


0.61 


1.00(5)(6) 


3-16 


0.27 


6.76 


0.007/0.05 


3 


0.58(4)(3) 


-10(2) 


0.56 


0.97(3)(3) 


3-16 


0.28 


6.76 


0.005/0.05 


3 


0.60(6)(4) 


-24(6) 


0.53 


0.98(3)(4) 


3-21 


0.56 


7.18 


0.031 


2 






0.64 


0.71(1)(4) 


5-25 


0.83 


7.11 


0.0124/0.031 


3 


0.47(6) 


-7(3) 


0.48 


0.64(2)(3) 


5-30 


0.49 


7.09 


0.0062/0.031 


3 


0.43(2) 


-22(3) 


0.40 


0.69(2)(3) 


4-30 


0.64 



TABLE XII: Results of two- and three-state fits to kaon propagators. 



It is worth pointing out that we fit these excited state masses in wall source propagators 
that were designed specifically to minimize the contribution of excited states. It is likely 
that analysis with other quark sources would further enhance our ability to resolve excited 
states. 

We note that the consistency of the excited K and it states with experiment indicates 
that there is no unphysical scale in these channels of length >2 lattice spacings. This is 
encouraging, since non-localities that might be introduced by taking the fourth root of the 
staggered determinant could show up here. 
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FIG. 11: Two-state fits to three-flavor kaon pseudoscalar propagators as a function of minimum 
distance included in the fit from the run with 10/g 2 = 7.09 and amug = 0.0062/0.031. The size of 
the symbols is proportional to the confidence level of the fit. Octagons and squares represent the 
two 0~ states, although, as discussed in the text, all of the confidence levels for this fit are so low 
that these symbols are extremely small. Standard size crosses are used for points where both the 
error bar and the confidence level are too small to be visible otherwise. These fits used -D m ax = 30. 

VII. CONCLUSIONS 



In this project we have calculated hadron masses including the effects of three flavors of 
dynamical quarks, using light quark masses down to 0.1 m s and lattice spacings of about 0.12 
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FIG. 12: Three-state fits to three-flavor kaon pseudoscalar propagators as a function of minimum 
distance included in the fit from the run with 10/g 2 = 7.09 and amug = 0.0062/0.031. The size of 
the symbols is proportional to the confidence level of the fit. Octagons and squares represent the 
two states; diamonds represent the oscillating 0~^~ state. -D m ax = 

30. 



and 0.09 fm. These quark masses are light enough that we are beginning to "see hadronic 
decays" in the sense that the lowest energy states for some quantum numbers may be two- 
meson states instead of a single particle. To the extent that we can reasonably expect, our 
spectrum results are consistent with the experimental hadron spectrum. One quantity that 
is sensitive to the effects of sea quarks is "J" , which is roughly the derivative of the vector 
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FIG. 13: Quenched kaon mass fit plot showing ground state and first excited state, with 10/<? = 
8.40, ami js = 0.0062,0.031, and L> max = 17. Fit to Eq. (J7J) without an oscillating state. 

meson mass with respect to the squared pseudoscalar mass[2^|. In particular, we plot 

J ~ 2{Ml-Ml) ■ m 
This quantity is plotted in Fig. HHJ which updates results from jlj , and also includes recent 
points from the CP-PACS/JLQCD collaboration^- 

Comparison of lattice results with the physical spectrum still requires extrapolations to 
zero lattice spacing and to the physical quark masses. In principle, the extrapolation to zero 
lattice spacing is straightforward — we expect errors proportional to a 2 g 2 . Extrapolation to 
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FIG. 14: Summary of fits of kaon propagators. Ground state and excited state kaon masses are 
interpolated to the correct strange quark mass. + parity partner state and n + K masses are 
uncorrected for comparison. For the Ks, ir and + entries in the legend, the symbol on the left 
represents fine lattice results, and that on the right the coarse lattice results. 

the physical light quark mass is more difficult. First, most of the hadrons decay strongly, and 
as we have seen for the ++ , and the + for nondegenerate quarks, simulations with light sea 
quark masses show the couplings to the decay channels. For stable hadrons the extrapolation 
to physical light quark mass involves chiral logarithms. Because of the remaining breaking 
of taste symmetry, fitting to the chiral logarithms requires that the continuum extrapolation 
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be done first, or simultaneously. 

In the case of the pseudoscalar masses and decay constants, taste violations have been 
included in the chiral perturbation theory, which makes possible a simultaneous extrapola- 



tion in lattice spacing and quark masses 



The small statistical errors on pseudoscalar 



masses and decay constants make this rather involved analysis necessary, but also make it 
possible. Work towards comparable extrapolations for some other quantities, such as the 
nucleon mass, is in progress. 

In the meantime, it is interesting to use a less sophisticated extrapolation to see how 
these lattice results compare with the real world. Figure ITH1 shows such a comparison, using 
a linear or quadratic extrapolation in the light quark mass and linear extrapolation in the 
squared lattice spacing. Since the difference between the strange quark mass used in our 
simulations and the correct value is roughly twice as large in the coarse runs as in the 
fine runs, the extrapolation in lattice spacing also largely corrects for the too-large strange 
quark mass used in the runs. (It is not entirely an accident that the continuum extrapolation 
largely takes care of adjusting the strange quark mass, since one of the largest reasons for 
the error in adjusting the strange quark mass was the neglect of order a 2 corrections in 
tuning the strange quark mass.) Note that the lattice nucleon mass plotted here is the 
linear extrapolation shown in Fig. HJ a proper chiral extrapolation is expected to lower this 
value. 

The spectrum results from these simulations with three dynamical light flavors are encour- 
aging. Clearly, however, considerably more work is needed, in particular on chiral extrapo- 
lations, before we can be confident that the calculations can produce accurate and precise 
results in all the channels that we have examined. Runs are continuing for m u ^ = 0.1m s on 
both coarse and fine lattices. 
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